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Abstract 

We consider the ground state of a trapped Bose-Einstein condensate in two di- 
mensions. In the mean-field approximation, the ground state density profile satisfies 
the Gross-Pitaevskii equation. We compute the leading quantum corrections to the 
density profile to second order in an expansion around the Thomas-Fermi limit. By 
summing the ladder diagrams, we are generalizing Schick's result for the ground state 
energy of a homogeneouns Bose gas to the case of a trapped Bose gas. 
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I. INTRODUCTION 



The remarkable achievement of Bose-Einstein condensation (BEC) of alkali atoms in 
harmonic traps B]-P] has created an enormous interest in the properties of dilute Bose 
gases (see e.g. Ref. |4[] and references therein). One of the basic questions is the condensate 
density profile as a function of temperature. At zero temperature and in the mean-field 
approximation, the condensate density satisfies the Gross-Pit aevskii equation. Until recently, 
the condensates were so dilute that mean-field theory gives a satisfactory description of the 
experiments. However, there are corrections from quantum fluctuations around the mean 
field and their relative importance grow like the gas parameter y/pd 6 , where a is the S- 
wave scattering length and p is the density. The sign and value of the scattering length 
are determined by the position of the least bound molecular bound states. The position of 
the bound states can be manipulated by applying external fields. In recent experiments, 
Cornish et al. were able to vary the scattering length a over a large range by applying a strong 
external magnetic field and exploiting the existence of a Feshbach resonance at B ~ 155 G. 
Values for \J pa 3 up to approximately 0.1 were obtained and should be sufficiently large to 
see deviations from the mean field in experiments. Hence it becomes important to be able 
to calculate the effects of quantum fluctuations around the mean field in a systematic way. 
Such an approach was developed in Ref. 0, where the effects of quantum fluctuations on 
the ground state of a Bose-Einstein condensate in three dimensions were calculated. In 
Ref. 0, the result was extended to arbitrary time- independent states, including vortices. 
The method is a combination of the Hartree-Fock approach and the Thomas- Fermi approach. 
The Hartree-Fock method includes all the leading quantum corrections to the mean-field 
equation for the density. The resulting equation is an integral equation, which is reduced 
to a local differential equation by applying a gradient expansion around the Thomas-Fermi 
limit. 

In a new paper, Blume and and Greene M, use a diffusion Monte Carlo method to 
calculate the ground state energy for different two-body potentials that generate the same 
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value for the scattering length a. It turns out that the ground state energy is independent of 
the actual form of the potential; it only depends on a. Moreover, it differs significantly from 
the ground state energy obtained from solving the Gross-Pitaevskii equation. Including the 
leading quantum corrections to zeroth order in a gradient expansion, the agreement with 
the result obtained by Blume and Greene || improves significantly. 

It is well known that a Bose-Einstein condensate in a two-dimensional homogeneous Bose 
gas only exists at T = p. At any finite temperature, the phase fluctuations destroy the 
condensate. This reflects the Mermin- Wagner theorem stating that there is no spontaneous 
breakdown of a continuous symmetry in a homogeneous system in two dimensions at finite 
temperature. 

The ground state energy density of a two-dimensional homogeneous Bose gas was first 
calculated by Schick |L0] . By summing up the ladder diagrams contributing to the chemical 
potential, he was able to obtain the leading correction to the mean field result. The result 
has later been derived by several authors jnUTjl . Very recently, a formal proof of the result 
by Schick was given by Lieb and Yngvason ]TB|jr^] . 

The Mermin- Wagner theorem does not apply to an inhomogeneous system and a 
condensate may exist. Bagnato and Kleppner jnj showed the possibility of BEC for 
T < T c « y/Nhu (where iV is the number of particles and u is the trap frequency) in 
a two-dimensional ideal gas confined by a harmonic potential. In this temperature range, 
the phase fluctuations on the scale of the size of a trapped gas is negligible, and there is a 
true condensate. In Refs. [0,0, it was shown that the possibility of BEC exists also when 
local two-body interactions are included. 

The ground state properties of a trapped Bose gas in two dimensions have been studied 



recently by a number of authors. In Refs. p!8|JT9|] , the Gross-Pitaevskii equation was solved 
numerically. However, the Gross-Pitaevskii equation is a mean-field equation and receives 



quantum corrections. In Refs. p^,I^,20|, it is argued that the correct equation for the 
density can be derived from an energy functional that includes the corrections obtained by 
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Schick |K| . 

However, since a trapped Bose gas is an inhomogeneous system, the Gross-Pitaevskii 
equation also receives corrections that involve gradients of the density. In the present paper, 
we calculate the leading quantum corrections to the Gross-Pitaevskii equation to second 
order in the gradient expansion around the Thomas- Fermi limit. We also apply the T- 
matrix approximation to sum up the ladder diagrams. The resulting equation extends the 
result in Refs. PIJS . 

The plan of the paper is as follows. In section II, we review the perturbative frame- 
work developed in Refs. || to calculate the leading quantum corrections to the the Gross- 
Pitaevskii equation. In section III, we calculate the self-consistent one-loop corrections to 
second order in the gradient expansion. In section IV, we briefly discuss the summation of 
the ladder diagrams and derive an equation that takes this summation into account. Finally, 
we summarize our results in section V. Calculational details are included in two appendices. 



II. PERTURBATIVE FRAMEWORK 

In this section, we discuss the perturbative framework developed in Ref. |J to calculate 
the leading effects on the ground state from quantum fluctuations around the mean field. 
The action is 

d h 2 



dt 



d 2 x V*(x,t) 



I d 2 x I rfV^*(x,t)^*(x',t)V (x-x')^(x,t)^(x',t) ) . 



(1) 



■0*(x,t) is complex field operator that creates a boson at the position x, V(x) is the 
trapping potential, Vo(x) is the two-body potential. In the following, we set h = 2m = 1. 
Factors of % and 2m can be reinserted using dimensional analysis. 

The action is symmetric under a phase transformation 



V>(x,t) -> e fa ^(x,t) 



(2) 



This ?7(l)-symmetry ensures that the density p and current density j satisfy the conti- 
nuity equation 

P + V-j = 0. (3) 



In the ground state, the current density j vanishes identically and the condensate has a 
constant phase. The i7(l)-symmetry can then be used to make the condensate real every- 
where. 

In the following, we approximate the two-body interaction Vo(x — x') by a delta function 
with strength g: 

^o(x-x') = g6(x-x!) . (4) 
We can then trivially integrate over x', and the action (|l|) simplifies to 

^-\g{r^f- (5) 

The quantum field theory defined by the action (|5]) has ultraviolet divergences that must 
be removed by renormalization of p and g. There is also an ultraviolet divergence in the 
expression for the density p. This divergence can be removed by adding a counterterm 
dp. Alternatively, one can eliminate the divergences associated with p and p by a normal- 
ordering prescription of the fields in Eq. (||). The coupling constant is renormalized in the 
usual way by replacing the bare coupling with the physical one. 

If we use a simple momentum cutoff M to cut off the ultraviolet divergences in the loop 
integrals, there will be terms proportional to M p , where p is a positive integer. There are 
also terms that are proportional to log M. The coefficients of the power divergences depend 
on the regularization method and are therefore artifacts of the regulator. On the other 
hand, the coefficients of log(M) are independent of the regulator and they therefore repre- 
sent real physics. In this paper, we use dimensional regularization to regulate both infrared 
and ultraviolet divergences. In dimensional regularization, one calculates the loop integrals 
in d = 2 — 2e dimensions for values of e where the integrals converge. One then analytically 
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dt / d 2 x ip* 



continues back to d = 2 dimensions. With dimensional regularization, an arbitrary renor- 
malization scale M is introduced. This scale can be identified with the simple momentum 
cutoff mentioned above. An advantage of dimensional regularization is that it automatically 
sets power divergences to zero, while logarithmic divergences show up as poles in e. In two 
dimensions, the one-loop counterterms for the chemical potential \x and the density p are 
quadratic ultraviolet divergences, while the one-loop counterterm for the coupling constant 
g is a logarithmic ultraviolet divergence. 

We next parameterize the quantum field ip in terms of a time-independent condensate 
v and a quantum fluctuating field ip\ 

ip = v + ip , (6) 

where the condensate v satisfies 

v = (V) • (7) 

Here and in the following, (A) denotes expectation value of the operator A in the ground 
state. Thus the expectation value of ip vanishes. The fluctuating field can be written in 
terms of two real fields: 

^ = -^(V>i+# 2 ) • (8) 
Substituting Eq. @ into Eq. (f|), the action can be decomposed into three terms 

S[lf)] = S[v] + Sfreefy'l, ^2] + S int [v, V>1, ^2] • (9) 

S[v] is the classical action 

S[v] = [ dt f d 2 x 



(li - V)v 2 - ^gv 4 + vV 2 v 



(10) 



while the free part of the action is 



Sfree[lpl,1p2] = dt d?X 



11 1 

(^^2 - ^1^2) + ^1 (V 2 - A 2 ) il) X + -V 2 V 2 ^ 2 



The interaction part of the action is 



(12) 



The sources in Eq. (|I2D are 



T 



fj, — V — gv 2 v + V 2 f , 



X = A + 



fi — V — 3gv^ 



Y 



fi-V - gv 2 



Z = -gv . 



(13) 
(14) 
(15) 
(16) 



Note that we have added and subtracted an arbitrary term |A ijj 1 in the action, effects are 
subtracted out at higher orders through the source X. By a judicious choice of A, one can 
ensure that X can be treated as a perturbation in the same way as the other sources and 
thus simplify calculations. We return to the choice of A in the next section. 

The propagator that corresponds to the free action Sfreet^ij ^2] i n Eq. (|TT| ) is 



D(u,k) 



uj 2 — e 2 {k) + ie 



k 2 —iuj 
ilu e 2 (k)/k 2 



J 



(17) 



Here k is the wavevector, u is the frequency, and e(k) is the dispersion relation for the 
Bogoliubov modes: 



e(Jfe) = kVW+A 2 



(18) 



The dispersion relation is gapless, which reflects the spontaneous breakdown of the £7(1)- 
symmetry (Goldstone's theorem). The dispersion is quadratic for large wavevectors and is 
that of a free nonrelativistic particle. The propagator is defined with an ie prescription in 
the usual way. The diagonal parts of the propagator are denoted by a solid and a dashed 
line, respectively. The off-diagonal parts are denoted by a solid-dashed and dashed-solid 
line, respectively. 



The quantum field theory defined by Eq. fl5|) has been decomposed into free and interact- 
ing parts. The quantum loop expansion is then an expansion in the dimensionless coupling 
constant g. The Hartree-Fock approximation includes all effects to leading order in g. 

The field equation for fa (and fa) is obtained by varying the action (0). Taking the 
expectation value of the equation for fa, we obtain the tadpole equation 

o = t + + l -z{^ - ^g(Mtf + tt)) ■ (19) 

The density can be written as 

p = (r*P) 

= v 2 + 1 -(rt) + l(i> 2 2 ), (20) 

To zeroth order in the loop expansion one neglects the expectation values in Eq. ( |i"9|) and 
the tadpole equation reduces to T = 0. We also neglect the expectation values in Eq. ( |20|) 
and we have p = v 2 . The tadpole equation then reduces to the Gross-Pitaevskii equation 
for the density profile: 

[^-V-gp^\+V 2 ^ = 0. (21) 

The last term in Eq. ( |l9|) only contributes at second order and beyond in the loop expansion. 
Thus to leading order in the quantum fluctuations, Eq. (|T9"D reduces to 

= T + h(tf)+ 1 -Z(^). (22) 

Eq. ( p2|) is referred to as the semiclassical tadpole equation ||. 

The sources X and Y depend on the condensate v. In order to obtain an equation for 
the density p, we invert Eq. fl20| ) so that we can eliminate v in favour of p in Eq. (p2|). 
The expectation values in Eq. ( pO"D are correction terms arising from quantum fluctuations 
around the mean field. These terms are suppressed by powers of g compared to the terms 
in the classical equations. Since we are only interested in the leading quantum effects, we 
can invert Eq. ( f20|) and expand it to first order in the expectation values: 
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(23) 



We can also derive expressions for the gradients of v by differentiating Eq. (\23\ 
1 



Vv = — 
Ap 

W 2 v ' 



4p 3/2 



4pV v ^+V^ 1 2 + ^ 2 2 )-v / pV^ 1 2 + ^ 2 2 )_ , 

V/ 2 v 2 Vp + \Vp v 2 v/p - 2 (v^) 2 ! + ^ 2 2 > 



(24) 



(25) 



Substituting the expressions ( |23| ) for t> and ( |25|) for V 2 v, the semiclassical tadpole equa- 
tion (E2h becomes 



= [li - V - gp] v/p+ V 2 v/p- g^l) + 



1 



2p 3 / 2 



-^~ v ~ 9P+ 7p v2 ^)^ + ^ 



(26) 



The last line in Eq. ( p6|) is proportional to the classical equation of motion. The corrections 
to the classical equation of motion are of order g and so are the expectation values of the 
quantum fields. Thus the last line is second order in quantum corrections and will therefore 
be dropped in the following. 



III. GRADIENT EXPANSION 

The expectation values in Eq. (^) are functionals of the sources X and Y. These 
functionals involve an arbitrary number of insertions of the sources in the loop diagrams for 
(ipf) and {ipl). They are non-local since the positions of the sources are integrated over. 
Thus the tadpole is an integral equation. 

Instead of solving the non-local integral equation (^), we would like to derive a local 
differential equation by expanding the expectation values (ipf) and (^f) in powers of the 
sources and their gradients. The Thomas-Fermi expansion is an expansion in powers of 
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£/R, where £ = 1/ y/gp is the local coherence length and R is the length scale for significant 
changes in the density. It is possible to expand the expectation values in Eq. (|2~2]) in powers of 
gradients of the sources, if the expectation values receive significant contributions only from 
modes with wavelengths that are of the order of the coherence length or less. This can be 
ensured by introducing an infrared cutoff that eliminates the contribution from wavelengths 
larger than the coherence length. We use dimensional regularization to guarantee that the 
effects of these modes are eliminated. If the results depend on the infrared cutoff, it indicates 
a sensitivity to length scales much larger than the coherence length ||. 

We can also expand the expectation values in powers of X and Y if the sources are 
at least either first order in the gradient expansion or in the coupling constant g. This is 
generally not correct, but it can be made true at a specific point xo by a clever choice of the 
parameter A ||. Since the expectation values of the quantum fields are of order g and we 
are interested in the leading quantum effects, we can use the classical equation of motion, 
T = 0, to simplify the sources: 

X(x ) = A 2 - 2gv 2 - ^V 2 v , (27) 
K(x ) = -U/ 2 v . (28) 

The source Y(xo) is already second order in the gradient expansion. We can also make 
X(x ) second order in the gradient expansion by the following choice || for the arbitrary 
parameter A 

A 2 = 2gp(x ) . (29) 

The difference between p and v 2 is higher order in g and can be neglected. Any other choice 
of A which is equal to Eq. (^) to leading order in g and in the gradient expansion, results 
in the same final equation for the density. 

The expectation values in Eq. fl26|) are now expanded in powers of X and Y and their 
gradients around the point x 
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(^ 2 ) = a + ai X + a 2 Y + a 3 V 2 X + a 4 X 2 + a 5 (VX) 2 + ... , (30) 
<^ 2 > = b + b x X + b 2 Y + b 3 V 2 X + 6 4 X 2 + 6 5 (VX) 2 + ... . (31) 

The expansions of (ipi) and (^2) include all rotationally invariant terms built from the 
sources and their derivatives. In Fig. p], the diagram that we need to calculate in order to 
obtain the coefficient qq appearing in Eq. (|30"D. This calculation is sketched in Appendix B. 




FIG. 1. Feynman diagram needed for the calculation of the coefficient in ag Eq. (|3[)| ) 



We will also need to calculate the gradients of {ip\) and (^f) in powers of the sources 
and their derivatives. Adding Eqs. (|30|) and ( |3~1~D and differentiating, we obtain 

V(^ 2 + V 2 2 ) = (a 1 + 6 1 )VX + ... , (32) 



V 2 <^ 2 + ^2) = (ai + &i)V 2 X + 2(a 4 + 64) (VX) 2 + .. 



(33) 



We next eliminate the dependence on v of the sources in favour of the density p. Since 
the sources only appear in the expectation values that are of order g, we can neglect the 
expectation values in Eqs. (p3|)-(p5|): 

u(xo) = v/p , (34) 
V^(x ) = Vv/p , (35) 
VMx ) = V 2 ^ . (36) 



The resulting expressions for the sources and their derivatives then become 

V 2 ^ 



X xo 



VX(x ) = -Ag^/p V^/p , 

V 2 X(x ) = -Ag ypV 2 yfp+{Vyfpf 
V 2 ^ 



(37) 
(38) 
(39) 
(40) 
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Using Eqs. (|30|)-(|33D for the expectation values and the expressions (p7|)-(|40D for the sources, 
the semiclassical tadpole equation reduces to 

= (p - V - gp) ■sfp + V 2 ^ - a g^/p 

+ 1 (oq + b + (2 ai + &i + a 2 )A 2 + 2a 3 A 4 ^ V v ^ 



2 v " ' p 

-\ (ao + b + ( 0l + h)A 2 - 2(a 3 - 2a 4 - 26 4 )A 4 + 4a 5 A 6 ) ^ ^} . (41) 

The coefficients and &j were calculated in Ref . . They can be expressed in terms of the 
integral J TOin , which is defined in Appendix A. We list them in Appendix B for convenience. 
The coefficients a and b are quadratically ultraviolet divergent, while a 1( a 2 , 6i, and b 2 
are logarithmically ultraviolet divergent. The remaining coefficients are ultraviolet finite. 
Note also that the infrared divergences in the individual coefficients cancel in the sum. The 
fact that the dependence on the infrared cutoff cancels, ensures that there is no sensitivity 
to length scales much larger than the coherence length to leading order in the gradient 
expansion. 

The counterterm needed to cancel the ultraviolet divergences in Eq. ([II]) is 



fc = £. (42) 



Using the values for the coefficients listed in Appendix B, we obtain 

9P 



i g 2 f 

P-V\^~p- gp^Tp - — 11 + log 



2M 2 

2 



PVP-^ 2 VP = 0. (43) 



Note in particular that the coefficient of (Vy 7 /)) in Eq. (|4i~|) vanishes. This is not the case 
in three dimensions || and we have no explanation for this cancellation. Eq. ( fi3|) has been 
derived at a specific point x . Since the point x is arbitrary, Eq. fl4~3|) must hold everywhere. 
The logarithmic term in Eq. (|43|) can be obtained from the one-loop result for the 



ground state energy density £ for a homogeneous Bose gas, obtained in Refs. [22,23[], by 
differentiating £ with respect to ^fp. The last term in Eq. ( fl3"D is a new result. 

This renormalization scale is completely arbitrary and physical quantities such as the 
density profile of the ground state must be independent of it. The requirement that physical 
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quantities be independent of M can be expressed in terms of renormalization group equations 
for the coupling constants in the Lagrangian. The coupling constant g in Eq. (|5]) satisfies 

s = m . (44) 

where the /^-function is a polynomial in g. Normally, these functions are known only up to 
a certain order in the quantum loop expansion. At the one-loop level /3(g) = g 2 : /47r. From 
Eq. (f44|), one can check that our result Eq. (|43|) is independent of the scale M. 

We have approximated the interaction by a delta-function which has zero range. In 
real systems, however, the interactions have a finite range a and so 1/a provides a natural 
ultraviolet cutoff. 



IV. T-MATRIX APPROXIMATION 



In this section we employ the ladder or T-matrix approximation, which takes an infinite 
number of loops into account. 

The diagrams we are summing and that give the leading correction to the mean-field 
term gp^fp in Eq. ( fi"3|) are shown in Fig. pj. The diagrams are not the conventional ladder 
diagrams since we use two real fields instead of a single complex field. The final result, 
however, is the same. The correction to the gradient term in Eq. (|43"D is obtained by 
summing the same diagrams with insertions of sources X, Y, and their derivatives in the 
leftmost loop in every diagram. 

The ladder diagrams are summed by solving an integral equation for the two-body 
T-matrix []IU],[23],|25|. In two dimensions, the summation amounts to replacing the bare inter- 
action g by an effective interaction 87rlog(2(ipa 2 ) _1 , which is density dependent [10,24,25 
In the dilute gas limit, this expansion parameter is small, and Eq. (43) is replaced by 



V 2 + p - V J-p - 8vr| log (2gpa l )\- l p^p - -| log (2gpa 2 ) \~ l = . 



(45) 



This is the main result of the present paper. The second term was obtained in Refs. OJT4IJ21] 
and follows from the corrections to the ground state energy of a homogeneous Bose gas 
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calculated by Schick ||10|| . The third term which is a correction to the gradient term is a new 
result. 

Up to corrections that are suppressed by powers of the small parameter log(2g pa 2 )^ 1 , 
Eq. ( f45|) can be derived from an energy functional E[$>], where |$|, is identified with ^fp: 



+47r|log(2^|$|V)|" 1 |$| 4 



l--|log(2^|$|V)|- 1 



(46) 



The energy functional Eq. fl4"6] ) generalizes the energy functional in Refs. P~5|j^D| ] to include 
a gradient term. 



V. SUMMARY 

In this paper, we have computed the leading quantum corrections to the Gross-Pitaevskii 
equation for a trapped Bose gas in two spatial dimensions. The method involves the trun- 
cation of two systematic expansions. The first is the quantum loop expansion which is an 
expansion in powers of the coupling constant g, and the second is a gradient expansion 
around the Thomas- Fermi limit. The result Eq. ( f43|) includes all leading order quantum 
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corrections to second order in the gradient expansion. 

The summation of the ladder diagrams changes the effective expansion parameter from 
g to 8tt log(2g pa 2 )^ 1 . The resulting equation includes a correction term that can be 
derived from the corrections to the ground state energy density of a homogeneous Bose gas. 
It also includes a new term which is a gradient correction to the mean field equation. 
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APPENDIX A: FORMULAS 



The loop integrals that appear in our calculations involve integrations over the energy u 
and the spatial momentum k. The energy integrals are evaluated using contour integration. 
The specific integrals needed are 

1 _ ?Y lVt+1 (-l)-l-3-(2n-3) 1 

2tt (u 2 - e 2 + ie) n 1 ' 2»(n-l)l e 2 ™" 1 ' 1 ' } 

*?_ = •/ nw+ i (-l)-l-3-(2n-3) 1 

2tt (lu 2 - e 2 + ie) n + l 1 ' 2"+%! e 2 ™- 1 ' 1 ' ' 

The momentum integrals are evaluated using dimensional regularization in d = 2 — 2e 
dimensions. Some of the integrals are infrared divergent or ultraviolet divergent or both. 
They can be written in terms of integral I m ,n, which is defined by 

/e 7 M 2 V r d d k k 2m 
m ' n ~ \ ' (2n) d k"(k 2 + A 2 )"/ 2 ' ( } 

Here, M is a renormalization scale that ensures that l m ,n has the canonical dimension also 
for d 2. 7 pa 0.5772 is the Euler-Mascheroni constant. With dimensional regularization, 
I mt n is given by the formula 
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(2vr) d V 4tt ; 2r(f) 
where tt d = 2it d l 2 /T[d/2] is the area of the ci-dimensional sphere. 
The integrals l m ,n satisfy the relations 

(d + 2m- n) J m> „ = nI m+2 , n+2 , (A. 5) 

A /m,n = Im—l,n—2 Im+l,n ■ (A-6) 

The first relation follows from integration by parts, while the second is simply an algebraic 
relation. 



APPENDIX B: COEFFICIENTS 

The coefficients needed to calculate (ipf) and (ipf) in Eqs. (|30| ) and (|3T| ) were calculated 
and listed in Ref. ||. For completeness, we list them below 



o-o = 2^1.1 > (B-l 

ai = ^2,3 , (B.2 

a 2 = — j!o,i , (B.3 

as = Wd [-10/5 ' 7 + 13/3 ' 5 ~ /l,3] ' (B ' 4 

3 

«4 = T^hfi , (B.5 
lb 

fl 5 = TTA [ 35/ 6,9 - 10J 4,7 + 3/2,5] , (B.6 

64a 

60 = , (B.7 

61 = -^0,1 , (B.8 

62 = ~/- 2 ,-i , (B.9 

63 = ^[2/3,6- 3/i,3 , (B.10 

64 = -^/i, 3 , (B.ll 
& s = [h,r+ 2/2,5 + /o, 3 ] • (B.12 
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The factor of d in the denominators in as, 05, 63, and 65 arises from averaging over angles 
in the momentum integrals. 

Finally, we outline how to calculate a one-loop diagram explicitly. For simplicitly, we 
compute the leading coefficient ao in the gradient expansion for the expectation value (ipf) 
in Eq. fl5Hp . The diagram is shown in Fig. [II and it reads 

/e 7 M 2 V fdur d d k ik 2 
a °~\^r) 7 2^7 (27r)^- e 2 + ie - (R13) 

After integrating over the energy u, using Eq. flA.l ), we obtain 

1 (e*M 2 Y r d d k k 2 

a 



2\ An J J (2ir) d e 
\h,i ■ (B.14) 
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